Digital computer control of polymerization process

ABSTRACT

Changes in component concentration, heat history, and other variables which cannot be directly measured in the successive reactors in a train of polymerization reactors, are followed by periodically numerically integrating, by a digital computer, for each reactor and for each of the variables, the equation WHERE X A PROCESS VARIABLE, E.G., CONCENTRATION, CONVERSION, ETC. I, AS THE FIRST SUBSCRIPT OF X, SIGNIFIES THAT THIS IS THE ITH OF I variables N, AS THE SECOND SUBSCRIPT OF X, MEANS THAT THIS IS THE VALUE OF THIS VARIABLE IN THE NTH ONE OF THE REACTORS ON THE LINE- THUS XI,1 MEANS THE VALUE OF XI IN THE FIRST REACTOR 64-1 ETC. F total volumetric flow rate V volume of the nth reactor T TIME DIFFERENCE UNDER THE REACTION CONDITIONS delta XI,N/ delta T OVERALL RATE OF GENERATION OR DEGENERATION OF XI,N UNDER THE CONDITIONS OBTAINING THE NTH REACTOR The resultant calculated values of these variables, together with directly measured values of other variables, are manipulated by the computer and used to adjust the rate of feed of reagents and other conditions of polymerization in the train.

[72] Inventor Robert C. Smith Primary Examiner-Malcolm A. Morrison Houston, Tex. Assistant Examiner- Felix D. Gruber [21] Appl. No. 48,954 An0rneysS. M. Clark and Willard L. G. Pollard [22] Filed June 22, 1970 [45] Patented Oct. 19, 1971 [73] Assignee The Firestone Tire & Rubber Company ABSTRACT: Changes in component concentration, heat his- Akron, Ohio tory, and other variables which cannot be directly measured in Continuation f application s the successive reactors in a train of polymerization reactors, 434 345 24 19 5 now abandoned are followed by periodically numerically integrating, by a digital computer, for each reactor and for each of the variables, the equation dx I )+8xi,n dz V z,n-l i,n at [54] DIGITAL COMPUTER CONTROL OF Where POLYMERIZATFON PROCESS at a process variable, e.g., concentration, conversion, 4 Clauns, 6 Drawing Figs. etc [52] 11.5. CI ..235/ 151.12, i, as the first subscript of x, signifies that this is the 208/159, 208/164, 260/767, 340/l72.5 ith of] variables [51] Int.CI G06f 15/46 n, as the second subscript of x, means that this is the [50] Field of Search 235/ 150.1, value of this variable in the nth one of the reactors 151.1, 151, 151.12; 260/1, 94.3-94.8,767; on the line-thus x 1 means the value of x,- in the 208/l59, 164 first reactor 64-1 etc. F total volumetric flow rate [56] References C'ted V= volume of the nth reactor UNITED STATES PATENTS 2 time difference under the reaction conditions 3,256,262 6/1966 Irvin 260/943 6m, 11 V 3,275,809 9/1966 Tolin et al" 235/15 1.12 overall rate ofgenerahon or degeneration ofx,, n under 3,296,097 1/1967 Lupfer 23s/151.12 x the ndltl0n5 Obtammg the reactor 3,321,280 5/1967 Trotter et al 235/15l.12X The resultant calculated values of these variables, together 3,351,430 11/1967 A h r 1, 235/151 12 X with directly measured values of other variables, are manipu- 3,378,483 4/ l968 Worrell et al 208/164 lated by the computer and used to adjust the rate of feed of re- 3,4l0,793 1 H1968 Stranahan et a1. 208/159 agents and other conditions of polymerization in the train.

I i l l ll To tit/M5025 /0 A22 Efifl A ZE K 4% AZZ A 4. I 20 "'l 111 111111111 if 2 u I lllll|1|1||||||1l|l|i 47 47 47 47 I 525C721: DIG/TAIL I 31 72: /ViUM/ZTIC SIGNAL tow/vzzs/o/v 72 1 II I l 20 ||Il|l|||1 mm 1 unnml nnllull g /4 COMM/T52 ll ill llllllllllll DIGITAL (lOll/lllPIJTER CONTROL OF POLYMERIZATION PROCESS RELATED APPLICATIONS This application is a continuation of my earlier copending application Ser. No. 434,845 filed Feb. 24, 1965, and now abandoned.

This invention relates to the polymerization of monomeric materials, and particularly conjugated diolefins such as butadiene and the like, and still more particularly to a method and plant for conducting such polymerization in an exactly controlled, reproducible manner.

Typically the conjugated diolefins such as butadiene are polymerized and copolymerized in continuously operating plants comprising a series of stirred reactors connected in cascade, each reactor continuously receiving discharge of the preceding reactor in the series and discharging continuously into the succeeding reactor in the series. The ingredients such as butadiene and styrene (hereinafter, for brevity the conjugated diolet'in monomer will be referred to as butadiene and the comonomer will be referred to as styrene, although it will be understood, as more fully explained below, that other monomers of these respective types may be used) emulsifying agent, catalyst, modifier, activator, etc. are fed into the first reactor in the series, the ingredients are progressively polymerized to form a latex as they pass through the successive reactors, and the resultant latex taken from the end of this series of reactors. A shortstop is added to arrest further polymerization, and residual monomers distilled out of the product. Control of this process is extremely difiicult, particularly in view of the long holdup time (on the order of 13 hours or so), and automatic control has not heretofore been successful. Accordingly the manipulation of the process has been left largely to the experience and discretion of human operators. In such a process the operator will have three fairly rigid targets in view. First there is the production schedule which is quite inflexible; so many tons of the particular type of polymer, and no more, have been ordered to come off the line for economic reasons, or there has been an obstruction in the coagulation, drying etc. equipment, and production must be curtailed. The operator generally takes care of this in a rigid manner by feeding the required amount of each raw material in an amount specified by the particular recipe to the system. Secondly there is the conversion to be controlled, i.e. the percentage of butadiene and/or other monomers which shall be converted to polymer during its passage through the line of reactors. This must also be rather rigidly controlled (1) so as not to overload the monomer stripping equipment and (2) more importantly, to maintain polymer properties such as Mooney viscosity, modulus, elongation and tensile strength. In order to achieve the desired conversion, the operator loads the ingredients into the system at relative rates which, in his judgment and bearing in mind the recipe and throughput schedules, will yield the desired results; but this is a very highly uncertain control due to the long holdup of the ingredients, and at best only approximate predetermination of conversion can be achieved in this manner. The final adjustment, which usually has to be rather extensive, is taken care of in a rough and ready manner by tapping off the product higher up, or lower down the line; e.g., if the product coming off at the 19th reactor is too high in conversion, the line is shortened by cutting out the 19th reactor and tapping off from the 18th or vice versa, if conversion at the 19th reactor is too low, the line is lengthened by cutting in an additional reactor after the 19th Much the same result can be obtained by injecting a polymerization arresting shortstopping agent at the reactor in the line where polymerization appears to have progressed to the desired conversion. Such control is obviously extremely cumbersome, since it involves physical cutting in and cutting out of major pieces of equipment on frequent short notice. This mode of operation also leaves greater or lesser proportions of the equipment idle at all times-It has also been proposed to raise or lower the temperature in the reactors in a train in response to deviation of the output from target values. As crudely carried out this is also very difficult and unsatisfactory, first because the long holdup makes such a control attempt subject to long oscillations, and also because it disturbs the properties of the product, which are not necessarily the same when the polymerization is carried out at different temperatures, even though the same final conversion is achieved. It will be noted that essentially the only effective control is maintained by adjusting conversion and temperature, and the effect of impurities, concentration of charge materials, system disturbances, pressure variations etc., are largely ignored.

Accordingly it is an object of this invention to provide a novel process and control for the polymerization and copolymerization of butadiene and other conjugated diolefms.

Another object is to provide such a process which may be readily controlled to produce polymeric products of uniform properties, particularly uniform and optimum conversion and Mooney viscosity, modulus, elongation and tensile strength.

Another object is to provide such a process which may be carried out with full and continuous utilization of all the reactors in the plant, without leaving a portion of the reactors unused as in the prior practice where control of conversion is achieved by cutting out a greater or lesser number of reactors at the latter part of the reactor series.

A still further object is to provide such a process in which the transition from one level of throughput to another may be readily effected without a decrease in product quality.

A further object is to effect the foregoing objects by means of controls imposed by a digital computer.

Other objects will be apparent from the discussion and description of the invention to follow. The invention will be described in connection with the accompanying drawings, wherein:

FIG. 1 is a flow chart of a plant for operation in accordance with this invention;

FIG. 2 is a diagram of a typical control for the flow of materials, responsive to and reporting to the control computer;

FIG. 3 is a diagram of a control analogous to that of FIG. 2 for the control of temperatures in reactors in the early part of the reactor train;

FIG. 4 is a diagram for a direct digital control of the temperatures in the later part of the reactor train; and

FIG. 5 is a graph showing the effects of proper and improper balance of catalyst and promoter in the reactor.

FIG. 6 is a graph showing the relation of conversion and Mooney MI. viscosity of a typical product at various modifier levels.

SYNOPSIS OF THE INVENTION Referring to FIG. 1, there is shown a synthetic rubber plant according to this invention. It is comprised of a controlled system of charge streams section 10 which are blended in a common line and charged to the train 12 of polymerization reactors whose temperatures are controlled by circulating methanol coolant through cooling coils inside the reactors. There is a digital computer 14 [in this specific instance, a TRW 330 computer, (trademark of the Thompson Ramo Woolridge Company) although other digital machines may be used] (a) which receives telemetered information through electric analog signal lines 17 and pneumatic signal lines 18 as to conditions which exist in the charge section and at various points of the reactor train 12; (b) which transmits telemetered control signals through lines 20 to the charge section 10 and to various points in the reactor train 12; and (c) into which the following information may be entered by the use of the operators console keyboard 16:

l. Management-dictated production requirements 2. Different type or types of production required 3. Laboratory analysis of current raw materials -4. Product properties, etc.

It will be understood that the measuring devices and control mechanism in the supply section and reactor train 12 may be of the pneumatic type or electronic type, whereas the computer 14 must operate upon digital values. In those cases where the measuring devices are of the pneumatic type there is accordingly provided a bank 19 of pneumatic-to-electrical analog converters receiving pneumatic signals from the pneumatic tubes 18 bringing information from the sensing devices in the equipment and converting these to analog electric current signals in corresponding electric lines 18. These electric analog signals in lines 13' are then converted, by a bank of electric-anaIog-to-electric-digital-signal converters 19', to digital signals in the corresponding line 18", through which they are transmitted to memory sectors in the computer corresponding to the variables whose values were measured and transmitted through the lines 18. These electrical to digital converters may or may not be contained within the computer.

Some of the telemetered measurement signals are in the form of electric analog signals coming through lines 17; these need only pass through the bank of electric analog-to-digital converters 19' to be impressed as digital signals through the lines 17' upon the memory of the computer.

Similarly, in those cases where the control mechanisms are of the pneumatic type, there is provided a bank 21 of electricdigital-to-pneumatic signal converters between the output transmitting lines 20' of the computer 14 and the pneumatic tubes 20 leading to the pneumatic type controller mechanisms in the chemical apparatus. It will be understood that the bank 21 of converters may carry out the conversion in two stages, first a conversion from digital to electrical analog output, and then from electrical analog to pneumatic output. The digitalto-electrical converters may or may not be contained in the computer. Of course if electronic measuring devices and electronic controlling mechanisms are used, no electrical-to-pneumatic or pneumatic-to-electrical converters are necessary.

It will be appreciated that to trace the individual lines 17, 18 and 20 from the computer to the actual points of measurement and control would tremendously confuse the drawing; instead therefore, the control and measuring lines in the supply and reactor trains l0 and 12 are traced for only short distances from their associated sensing devices and actuators; they are then terminated with one of the symbols or indicating that they are respectively continuations of the lines 17, 18 and 20 leading to the computer. The arrows on these symbols indicate the direction in which information or control signals are being transmitted.

Of course, if electronic measuring devices and actuators are used, no such conversion is necessary. In the drawing, the control and communication lines l7, 17, 18, 18', 18'', 20 and 20 are shown as dotted lines to distinguish them from material transport lines; arrows on these lines indicate the direction in which signals are transmitted, i.e., either to or from the computer 14. It is understood that there are many means of transmitting signals from the computer to the control mechanisms such as electrical analog, pneumatic or digital.

Before getting into the details of the operation of the plant, a word on the overall philosophy of the control of this invention will be found helpful. In general, the principle is to set up inside the computer 14 a mathematical model of the entire plant or of such section thereof as is involved in the control required at the moment. For the reactor train, the control model maintains in real time within the computer memory, in an array of memory sectors, hereinafter designated the model array the numerical values of all variables, which relate the current status of the material, in all the reactors which are in the reactor chain.

Thus at any time the computer memory may be interrogated to determine a value which represents the concentrations of any component or value of any other variable or parameter in any given reactor.

The values of these variables are obtained from analog inputs (such as flow rates and densitometer readings), from values entered through the operator's console (such as activator solution strength, target conversion desired and production rate to be maintained), and from values (parameters) calculated in the mathematical models.

At frequent intervals (say 5 minutes; the computer could perform this several times more frequently but this frequency is adequate) the computer performs the manipulations of the mathematical model which comprises the following general steps:

I. A simulation phase" updates the values in the model array" to the present conditions of all the variables in all of the reactors which are in service.

The values are either actually measured via lines 17 and 18; or entered through the operator's console since the model program went through its last cycle; or values calculated (as in the case of parameters) the last time the mathematical models were carried out.

2. A target rate of reaction is calculated which will achieve the desired conversion when the charge material at the charge end of the reactor chain exits the final available reactor in the reactor chain 12. The existing conditions of desired production rate and the number of reactors in the reactor chain which are considered available for service are considered in these calculations.

3. It calculates from the actual conversion measured by a radiation densitometer 96 at about the middle of the line, and adjusts to a standard temperature (say 0 F.) the average reaction rate prevailing in the reactor line up to the densitometer 96.

4. It then adjusts the numerical value of the activator rcquiremenLs in the memory section proportional to the discrepancy between the target and actual rates of reaction plus adjustments for relative changes in certain measured or calculated variables (such as soap pH for example).

5. It then performs similar calculations to adjust memory values of certain other variables (feed of supplemented catalyst at 92; temperature ofeach reactor in the reactor train past the midline of conversion determination) to provide a final trim to the target conversion.

6. Therealter the computer examines the memory sectors of all of the imposed and assumed operating values, and imposes these, as set points. through the lines 20, upon the mechanically operating control elements of the actual physical plant.

DESCRIPTION OF THE CHEMICAL PROCESS EQUIPMENT EXCLUSIVE OF CONTROL Initial Reactant Feed Section 10 This comprises a plurality of pumps 22, 24, 26, 28 and 30 drawing from sources of respectively butadiene, styrene, aqueous soap solution, activator and catalyst as indicated by the legends on the suction lines thereof and discharging through feed regulating control valves 32, 34, 36, 38 and 40. Likewise a source of mercaptan modifier is provided comprising a blowcase 42 pressured with nitrogen and containing the modifier which is discharged through the flow regulating control valve 44. The set points of the flow regulating control valves 32, 34, 36, 38, 40 and 44 are governed by telemetered pneumatic analog signals from the lines 48, 50, 52, 54, 56 and 58 indicated by the symbols Q? as extending back at 20 and 20 to the memory sectors in the computer 14 devoted respectively to set points for the feed of butadiene, styrene, aqueous soap solution, activator, catalyst and modifier. All of the flow regulator valves 32, 34, 36, 38 and 44 discharge into a header 60 leading to the intake 62 of the first reactor 64-1 in the reactor bank 12. The catalyst control valve 40 differs from the other valves in that, instead of discharging into the header 60 it feeds through a line 101 into the bottom of the reactor 64-2 in the reactor chain. Actual polymerization does not set in until the catalyst is added, and the first reactor 64-1 is devoted to precooling and mixing the ingredients. Alternatively, the line 101 could lead to any of the first three reactors, depending on the need for precooling. Thus at all times the various ingredients of butadiene, styrene, etc. will be fed into the headers 60 and reactor train 12, each at a rate in proportion to the numerical values currently in the corresponding memory sectors in the computer 14.

Details of Control Linkages At this point it will be of interest to digress and consider the details of the linkage between the computer and the materialsflow control actuators upon which it imposes set points. FIG. 2 shows details of the linkage involving the butadiene valve 32. The computer generates its commands in the form of digital impulses in the electrical line 102. These are transmitted to a digital-to-electric-analog converter 104, (one of the bank of such devices indicated at 21 in FIG. 1) which is a conventional device involving a servo adjusted potentiometer which changes its position whenever a signal differing from previous signals is transmitted to it, so as to generate a DC current in the line 106 proportional to the value of the digital signal. This in turn actuates an I/P cell 108, which is a conventional device which generates pneumatic pressure in the pneumatic line 43 proportional to the DC current in the line 1106. This pneumatic pressure is transmitted via the line 48 to a flow rate controller 106. In the butadiene supply line 110 there is provided a flow measuring orifice 112, the pressure drop across which activates a DP. Cell 114 which is a conventional device which generates a pneumatic pressure in the line 116 proportional to the square root of the butadiene flow in the line 110 and transmits this pneumatic pressure through lines 118, 119 to the flow regulator controller 100. This flow regulator controller compares the set point imposed by the line 48 with the flow indicated by the pressure in the line 113, 119, and if they are discrepant, changes the air pressure in a line 122 leading to the actuating diaphragm 124 of the valve 32, in such sense as to increase the butadiene flow if this is indicated to be deficient or to decrease this flow if excessive. Thus the flow of hutadiene in the line 110 should follow reasonably closely the value demanded by the original digital control signal generated at 102 by the computer 14.

However, the flow through the valve 32 may not follow the control signal exactly. Accordingly, in order that the computer may be informed of the actual flow, the flow signal line 114 branches to the return information line 49 which goes to an 1/? cell 126 which converts the pneumatic flow signal to a DC electric current in the line 120 proportional to the pneumatic signal. The current in the line 126 is converted by an analog to digital converter 129 into a digital electric signal in a line 1310 and transmitted to an appropriate sector of the memory 132' in the computer. It is this sector of the memory to which computer will refer for the value of the flow of butadiene in conducting calculations as hereinafter described. Calculations conducted by the computer will thus all be based upon actual values of flow, etc. as reported by the sensing device, rather than upon the values of set points imposed.

It will be understood that the same general system is applied to all of the control actuators in the system, for control materials flow; i.e., for each of the valves 34, 36, 38, 40, 44 (and also 78 not yet mentioned) there will be respectively pneumatic control lines 50, 52, 54, 56, 58 and 79 similar to line 48 for valve 32; back-reporting lines 51, 53, 55, 57, 59 and 81 similar to line 49 for valve 32; and associated actuating and measuring equipment similar to that shown for valve 32 in FIG. 2. The Reactor Bank 12 This comprises a series of reactors 64-1, 64-2, -3, -4 and so on down to 64-21, of glass-lined steel or other corrosion-resistant construction. Each of the reactors is provided with a rotary stirrer 66 and a cooling coil 68 to which cold methanol is admitted through a temperature regulating control valve 70 which is governed by the temperature reigning inside of the reactor, so as to admit cold methanol whenever the temperature in the reactor rises above the set point for the instrument. These auxiliaries are shown sketchily in FIG. 1 for reactor 64-2, and more in detail in FIG. 3 and 4 discussed below; it will be understood that each of the reactors 64-1 to 64-21 will have a temperature control of this general sort. The set point for each of the instruments is received from the computer, and the actual temperature of the reactor contents transmitted back to the computer as analog input. The temperature controls on the reactors up to the midline point are controlled by a common pneumatic analog control line 46, i.e., it is attempted to keep the temperatures in these reactors all at the same level; and the true temperature of the contents of each individual reactor is'transmitted back to the computer via a resistance thermal detector signal line 47 as described in connection with FIG. 3 hereinbelow. Reactors 64-11 to 64-21 each have an individual direct digital control for temperature as described in connection with FIG. 4. By direct digital control it is understood the computer sends an electrical signal which is transduced to pneumatic and applied directly to the control valve, thus bypassing the conventional local controller. It will also be noted in FIG. 4 that an independent temperature system is provided which records individual reactor temperatures and is useful in controlling individual reactors in emergency circumstances when the computer is not in control.

Reverting to the flow of materials, each reactor in the series is connected at its top to a pipe 72 leading to the bottom of the next reactor in the series; thus the ingredients introduced into reactor 64-1 at 62 are directed to pass successively through the reactors 64-1, 64-2, 64-3 and so on and ultimately down to reactor 64-21. The reactants polymerize progressively during their passage together through the reactors, the end product coming from the discharge pipe 74 at the top of the last reactor 64-21 being a latex of a copolymer of butadiene and styrene containing greater or less proportions of unreacted monomers. Typically at this point it will be desired to have about 60 percent of the butadiene reacted to form a copolymer and the remaining portion of 40 percent still unpolymerized; hereafter this will be regarded as the target conversion desired, although it will be understood that greater or lesser conversions could be set up as targets. The latex passes through the pipe 74 to a junction 76 at which a shortstopping agent is injected through a line 78 to stop all further polymerization and the killed latex is then passed through a line 80 to the monomer flash-off tanks and finishing equipment of the plant, none of these latter items being shown. The rate of flow of the shortstop is under control of the computer to the extent necessary to kill the reaction.

Radiation densitometers 96 and 94 are provided in the line 72 just ahead of reactors 64-11 and 64-13 respectively, and transmit this information to appropriate memory sectors in the computer 14. Two densitometers have been provided to impart flexibility to the operation. Normally densitometer 96 is in service and densitometer 94 idle; however, if is desired to perform maintenance on densitometer 96, densitometer 94 may be used in its place.

The last few reactors 64-11 to 64-21 have individual temperature controls 47. The temperatures are controlled individually, based upon mathematical models in the computer such that the latex in each individual reactor will achieve target conversion when it exits the final reactor at 80. This temperature control serves to correct any errors that may have occurred in the previous conversion control models. The temperature control consists of three parts: first a feed-forward portion based upon information obtained at the midline densitometer 96 that determines the temperature requirements of the latex currently at the midline point; secondly a compensation or correction applied to each reactor that corrects for possible poor response of the reactor temperatures to a temperature set point change; and thirdly a feedback correction based upon the error between the conversion measured at the end-of-line densitometer 102 and the target conversion. The feedback portion of temperature control in effect corrects for errors which may have occurred in the feed-forward portion of the temperature control.

A metering pump 92 is provided for injecting additional catalyst at a point between the midpoint of the train and the end thereof to reactor 64-14. The stroke of this pump is under the control of the computer 14 via a control line 86. The purpose of the metering pump 92 is to inject additional catalyst if the computer mathematical models determine that the reaction will not achieve the target conversion due to imbalance of the ratio of catalyst and activator notwithstanding the theoretically correct charge injected by the regulated feed section 10.

Referring now to details of the temperature controls for the reactors, and first to FIG. 3 showing reactor 64-2 which is typical of the early reactors in the line, namely reactors 64-1 to 64-10; as shown, the control comprises the pneumatic line 46 common to all of said early reactors and carrying a temperature set point pressure signal generated by the computer, which line 46 is connected by a branch 119 to a temperature recorder and controller 124. A temperature measuring bulb 120 is immersed in the contents of the reactor, and transmits the temperature in the reactor through the capillary tube 122 to the temperature recorder and controller 124. This temperature recorder and controller compares the temperature indicated via the capillary 122 with the pneumatic set point pressure received via the branch line 119, and if the temperature indicated via the capillary 122 is greater than indicated by the set point, the controller 124 alters the air pressure in a line 126 going to the diaphragm 128 on the valve 70 in such sense as to increase the flow of cold methanol through the valve 70 into the cooling coil 68 immersed in the contents of the reactor thereby lowering the temperature of the reactor contents. Conversely, if the temperature indicated via the capillary 122 is lower than the set point temperature, the controller actuates the valve 70 to decrease the flow of cold methanol, whereby the heat of reaction is pennitted to raise the temperature of the reactor contents. In this way, the temperature of the reactor is kept closely equivalent to the value present in the corresponding memory sector of the computer 14.

The actual temperature in the reactor will be subject to excursions, and in order that the true temperature may be at all times communicated to the computer, there is provided a resistance thermometer element 130 inside the reactor which generates an electrical analog signal in the line 47 which is transmitted as indicated at 17 to the appropriate memory sec tor of the computer 14.

The temperatures in each of the later reactors in the train, viz. 64-11 to 64-21 are individually controlled by a semidirect digital control as shown in FIG. 4, which is typical for all of these reactors. In FIG. 4 is shown the reactor 64-16 having a rotary agitator 66, cold methanol valve 70 and cooling coil 68 similar to the arrangement in FIG. 3. The temperature prevailing in the reactor is measured by a resistance thermometer 140, transmitted as an electric analog signal through the line 47 to an analog/digital converter 144 which transmits corresponding digital signals to a memory sector (a) in the computer. There is in the computer a temperature set-point memory sector (b) containing the temperature desired in the reactor and also a further memory sector (c) which sends out (via the line 148, digital-to-electric-analog converter 150 and electric-analog-to pneumatic pressure converter 152, through the line 154) to the diaphragm 128 of the valve 70, a pneumatic pressure corresponding to the number in the memory sector (c). The computer periodically compares the memory sectors (a) and (b). If the temperature in the reactor as recorded in memory sector (a) is higher than the desired value in memory sector (b), the computer changes the memory sector (c) in such sense as to open the valve 70 wider and to cool off the batch. Conversely, if the temperature in the reactor is too low, the valve 70 will be cut back so as to allow the heat of reaction to increase the temperature in the reactor 64-16. Thus the reactor temperature will be at all times kept approximately equivalent to the numerical value in the memory sector (b), and a continuously up-dated value of the actual temperature in the reactor will always be present in the memory sector (a). It will be understood that each of the reactors 64-11 to 64-21 will be provided with its own independent system similar to that described above; each reactor will be temperature controlled in accordance with its own respective individual temperatures set-point memory sector independently of all of the other reactors.

THE MATHEMATICS AND PROGRAMMING OF THE COMPUTER 14 where x= a process variable, e.g., concentration, conversion, etc. i, as the first subscript of x, signifies that this is the ith of 1 variables n, as the second subscript of x, means that this is the value of this variable in the nth one of the reactors on the linethus x means the value of x, in the first reactor 64-1 etc. F total volumetric flow rate V= volume of the nth reactor z= time difference under the reaction conditions (fir /8!) overall rate of generation or degeneration of x under the conditions obtaining in the nth reactor In equation (1), the first item on the right-hand side will be seen to be purely the result of materials balance. The second item on the right represents the generation or degeneration of the variable x (e.g., the increase in conversion, or some mathematical representation of the thermal history of the material in the nth reactor such as HDT," a concept which will appear below; this second term may be a more or less complicated function of several variables-temperature, residence time in the reactor, conversion, etc. In the simulation subroutine, for convenience, since the form thereof is identical for all variables, the computer first runs down all the variables and corrects their values solely for materials balance since the last simulation, by a routine which, in terms of Fortran, would go as follows, assuming (as will be seen below in Table II) that 24 variables are to be tracked, that we have all 21 of our reactors on line and that the time since the last updating and simulation was 5 minutes:

TABLE I DOlOl=Bl, 24

DOlON=Bl, 21

10 X(I,N)=X(I,N)-l(F/V)" (X1,N-l )-X(I,N))*5.0 The notation is the same as for equation (1 except that capital letters must be used in place of lower case letters. This subroutine evaluates the first term of equation (1) 24 2l times-in parallel for each of the 24 variables and in series for each of 21 reactors. In the claims hereafter the sequence of operations will be expressed in terms of the FORTRAN programs of Tables I and 11. However, it will be appreciated that the updating cycle interval of 5 minutes is merely chosen for convenience and that other time intervals could be used; accordingly, for generality, the time interval factor 5.0" in these equations will be expressed as DTlME" in the claims. Likewise the assumption of 21 reactors are on line is purely arbitrary and accordingly the index 21 for N on the DO loop is replaced by the symbol NUMBER" in the claims. Also, for simplicity, the claims will consider only the modification of one variable at a time, say the concentration, so that the double-subscripted variable X(I,N) need only be a singlesubscripted variable X(N); and the first "DO statement of Table I may be omitted. So condensed, the program of Table I appears in the claims as follows:

TABLE I-A D l0N= 1,NUMBER l0 X(N) X(N)+(F/V)*)X(Nl )X(N))* DTIME After execution of the program of Table I, the simulation subroutine then passes on to those variables for which the second term of equation (1) is not zero. For instance, there has been developed in connection with this invention a thermal history figure HDT which is associated with the latex stream at any point, such that:

where =2 l8, a factor serving the purpose of temperature adjustment roughly in accordance with the Arrhenius equation T= temperature in the reactor in degrees Fahrenheit t= time of subjection of the reactants to temperature T This quantity l-IDT may be considered as the total reaction time which the latex has undergone, referenced to 0 F. This quantity is blended forward with the flow in accordance with the subroutine of Table I above. It is then further updated for generation as well as mixing by a subroutine represented as follows:

TABLE II DO ll N 1,21

ll HDT(N) HDT (N) H(N) 5.0 The 5" in statement 11 represents the 5 minutes elapsing since the last simulation cycle updating. The various T(N)'s in the second statement were of course present in memory by report from the temperature sensing elements 130 in the field via the lines 17. Similar further corrections are made in all other variables for which, in equation (2) above, the second term on the right-hand side is not zero, this being in accordance with common methods of numerical integration. Certain of the claims will be directed to the updating of HDT alone, without consideration of other variables. In such case HDT(N) is a specific instance of X(N) of Table I-A, and for the purpose of the claims I-IDT in Table II is replaced by X(N and again the time interval of 5.0 minutes is replaced by DTIME"; and the number of reactors is expressed as NUM- BER instead of 21. So condensed, the program of Table II is as follows:

TABLE ll-A D0 ll l,NUMBER H(N)=2 (T( ll X(N) X(N) H(N) DTIME As stated above HDT is updated by first updating for materials balance by executing Table I (simplified as I-A) to update l-IDT for materials balance and then executing Table II (simplified as Table II-A) to update HDT for generation in the reactors.

An exemplary preferred machine-language program for updating the variables for materials balance in accordance with Table I is given hereafter in Table IV, and a similar machinelanguage program corresponding to Table II is given hereinafter in Table V.

The above is the general type of routine used for the summations and integrations to be discussed below, and, for brevity, the methods of programming will be taken for granted and the process discussed in terms of mathematics apart from programming.

A tabular summary of the channels-available in the computer M for simulation and tracking, and of the variables simulated and tracked, is given herewith in Table III.

TABLE III I. Butadiene impurity 0 l 2. Butadiene impurity 02 3. Butadiene impurity 03 4. Butadiene impurity 04 5. Modifier 6. Activator 7. Soap (parts) 8. Soap pH 9. Butadiene impurity 06 I0. Styrene impurity 02 Percent oxygen in butadiene vapor Reactor rate Density of stream entry reactor 64-1 Reducing power of activator solution Styrene impurity 03 Calculated conversion DT (reaction time up to reactor in question) HDT-temperature history referenced to 0 F. Catalyst Unused inch inch inch inch inch inch 27. Midline conversion of latex post reactor 64-10 28. l-IDT* 29. Midline catalyst blended to end of line 30. Incremental catalyst parts from Some of these are used in the computer routines described below, and others are simply recorded for analysis of the system and possible inclusion in future, more refined, calculations.

Having completed the simulation subroutine, the computer then goes on to decide on any adjustments to be made, as follows:

A, Control of Activator Flow at 38 l. The temperature in all the first reactors 64-1 to 64-13 is assumed to be 42 F.

II. The temperature in the later remaining reactors is assumed to be 44 F.

Holdup time per reactor (1-) is at target production rate levels.

Assume initially, and adjust on each successive Sminute cycle, a reaction rate until the conversion rojected according to the steady-state conversion equation. is equal to target conversion, To this end the computer first evaluates the equation.

RTnew iold the numerical value 60" in equation (5) is usual value of the management-dictated percentage conversion; this could, of course, be some different figure. The R soon converges on the correct rate of reaction which the latex at the first actively polymerizing reactor 64-2 will require to maintain target conditions when it arrives at the end of the line.

rcactur 64-10 Hdt (6) I'm, CALI.

reactor -2 wherein r apparent reaction rate up to midline, referenced to C is the actual conversion, at midline, measured by the densitometer 96, and

the integral expression is that of equation (3).

The machine then determines what activator level is required to maintain the target rate of reaction (R per convergence of equations (4) and (5) above) by the equation. 7. Ao=AML*(R /r wherein A0 level of activator required to come from the feed line 60 AML activator level at midline, previously determined in the simulat on routine From A0 and the throughput, the computer calculates the rate of flow of activator required through the valve 38, and inserts this value in the memory sector which governs the set point of valve 38. This completes the adjustment of activator supply, which is the initial and principal control maintained by the system of this invention.

While the invention has been described as controlling the reaction by the feed of activator, it will be appreciated that the control could also be exercised for other variables, e.g., soap concentration, catalyst feed concentration, etc.

B. Control of Temperatures in the Past Midline Reactors At the midline densitometer 96 the first opportunity arrives to determine the conversion and calculate the rate of reaction. Using this information, the computer will calculate the temperature that the latex at the effluent of the reactor 64-10 will require to maintain target conversion when this latex exits at 80 many hours later. 5 minutes previously (for example) the latex may have required a different temperature due to a different rate of reaction. This is handled through the materialsbalance portion of the simulation subroutine. This enables us to track the temperature requirements in real time, and provides, by interrogating memory, the temperature requirements in each reactor from midline to the end-of-line reactor.

To determine the temperature requirements of the latex at reactor 64-11 and later reactors the computer employs a trialand-error solution; first assuming a temperature and using the steady-state form of the conversion equation (4) to determine what the conversion will be on exiting at 80.

l. The computer assumes a value H, for the reactor 64-11,

then evaluates the equation.

reactor 64-11 wherein r is the midline value of reaction, determined per equation (6) above referenced to 0 F. 1- holdup time per vessel in hours I]. The computer then corrects H, for the next trial (5 minutes hence) by the equation.

where C is the measured conversion calculated from the reading of the midline densitometer 96. 60" is the desired target conversion, which could, of course, be some other arbitrary figure. Equations (8) and (9) quickly converge after a few trials on a value of H, that will give the target conversion 60.

Ill. The computer then determines the thermal treatment HDT" of equation 10. HDT ,,,,,,,,=H ,*NO*-r where N0 number of downstream reactors and this is the required future reaction time referenced to 0 F., that the latex in reactor 64-10 must undergo before exiting at in order to obtain 60 percent conversion.

IV. The computer then determines the total, thermal history referenced to 0 F., that the latex must have on exiting at 80, by evaluating the equation. 11. HDT"( l0)=HDT ,,,,,,,,+HDT(l0) wherein HDT*(10) total reaction time, referenced to 0 F., that the latex must have on exiting at 80 in order to hit target conversion V. The simulator routine blends HDT* associated with the reactor 64-10 downstream into the successive downstream reactors by a subroutine similar to that of Table I, but with the DO loop indexed from 11 to 21, operating on the equation. 12. HDT*(N)=HDT(N)+R(HDT*(N-l )HDT*(N)) where R Flow At/(volume of reactor N+Flow*Al) volume fraction blended V1. The computer then calculates the temperature requirements of each reactor by the equation.

wherein HDT*(N) reaction time referenced to 0 F. required for the latex in vessel N reactor N-l HDT(N 1) reactor 2 Hdt=t0tal reaction time, referenced to 0 F., that the latex in reactor N1 has been subject to N0 number of vessels remaining in which to trim temperatures r= holdup time per vessel, hrs.

a feed-back connection for the fact that the feed-forward model is in error DHDT= wherein TS (N) is the set-point temperature for the Nth reactor, and

is inserted in the memory sector that controls Nth temperature of the Nth reactor. This completes the temperature trim portion of the control. C. Control of Catalyst Addition Referring to FIG. 5, the time-conversion polymerization curve has a different shape under different activator-catalyst ratio relationships. The central curve (2) represents the conversion equations (4) etc. upon which the control of this conversion is based. The simulator routine maintains the calculated conversion in each vessel in memory in real time by the equation. l5. C(N)=C(N)+R(C(N-l)-C(N) (Flow Portion) l6. C(N)=C(N)+r(N)*H(N)*X(N)* 5. (Reaction Portion) manipulated per the simulation routines of Tables I and II above. Thus assuming that the model-calculated conversion at the end of the line is at the point C and the actual conversion as measured by the densitometer is at A, we know that this represents conditions of deficient catalyst and conversely if the measured conversion is at B, we know that catalyst is excessive. If the measured conversion falls on point C, we know conditions are correct.

Bearing this in mind, the catalyst charged at the beginning of the line at 101 is handled on a feedback basis. We have one constraint, namely that the minimum catalyst level for a given activator level should preferably be:

17. CAT ,,,=(9.0*A+20.0)A where CA T minimum molarity of the catalyst A molarity of the activator At the end of the line:

( i l 2 1 of m (EJJ CM (EaL ratio of requirements to minimum (see eq. 17) then CAT() =X*(9.0 *A (2) 20) *A (2) required parts of catalyst to be charged at reactor 64-2 where CAT(0) =catalyst required to be charged at 4X c ato lsyg 'jnjststqt'klm. The computer then changes the set point of the catalyst valve 40 in the indicated ratio. Incremental Catalyst Added at 92 The amount XCAT of catalyst required to force CC (EOL)=CM(EOL) is calculated by the machine sequenc- (a) is negative no XCAT is required. (b) is positive XCAT is required. NXC AT refers to the reactor where incremental catalyst is to be added. 5 X (0) refers to the charge at 101.

Thus an error in the CC-CM not only corrects the charge catalyst at 101, but also initiates corrective action by means of incremental catalyst at 92.

A second portion of the incremental catalyst program in- 1 volves prediction of whether it is possible by temperature trim alone to control the latex at reactor 64-11 to maintain conversion at target when the latex exits at 80. If the temperature is insufficient, within permissible limits (say 40-48 F.) incremental catalyst is added by increasing the stroke set point for the metering pump 92. THe equations used are The net of the incremental catalyst subroutine is to calculate the (XCAT) TOTAL to 1) correct both for improper balance of catalyst and (2) make good the inadequacy of feed-forward temperature to maintain target conversion, the sum being XCAT(TOTAL)=XCAT+ AXCAT the memory sector containing X CAT (TOTAL) is connected through the line 86 to set the stroke of the metering pump proportional to this figure.

D. Miscellaneous Control Loops Utilizing simulation and target adjustment principles as described in connection with the preceding control loops, the computer makes use of the following equations:

1. Shortstop control at 76.

The shortstop addition is controlled by the valve 78, whose set point comes from the computer via line 79. The shortstop is made proportional to the simulated amount of catalyst by evaluating in sequence the equations:

26. SHORTSTOP (EOL) CATALYST (EOL)"G 27. Check against a predetermined set of control limits 28. Calculate R value ratio of calculated parts in Equation ing through the equations. 2 to target parts XML): CATiElm/ACT (20) This R value is used to adjust the set point of the valve 78. gig (506M) G 2. Modifier The laboratory-determined Mooney MS, value of the (AcTME'Lfi-m) product and the corresponding line conversion are entered 64r2/XLL (EH1) periodically (say every 2 hours) via the console keyboard 16. XWXCATF CATWXCAT) From actual plant experience data a correlation of Mooney 64-14 64'l4/ACT(NXCAT) MS, vs. conversion at different levels of modification was 64-14 65 determined. FIG. 6 shows the general form of the correlation, LLtzvxcnF) the approximately straight lines being contours for various a 64-14 parts of modifier. From the slope of the lines the relationship X (NXCAT) 4 %Z I VXCAT) between MS-4 and conversion is established. In the sequence [XO6@;,1=/X0(NXAD] *XLMNXCAHMCTWXCAT) of equations below, first a conversion is calculated to cor- Q 41 5PM 4 '1'4 21 respond to the keyboard-entered laboratory MS-4. Next this conversion is compared with the actual conversion that exits t ;m in the process. From this an error is established. A portion of 64' 2 the error is combined with the error that existed the last time the MS-4 conversion values were entered. This blended error is used to adjust the present charge of modifier (M in equation 32). The sequence of processing the equations is as follows:

29. Cale. conv. (Lab. MS-4 K )"G,+K

30. E (error) (Lab. conv. calc. conv.)

33. Modifier parts present modifier charge M,

34. Check modifier parts against control limits 35. Calculate: R calculated parts (in (33) )/target parts 36. This R is used as basis for the calculations of the set points for output to the control mechanism 3. Production Rate Changes Each type rubber has its own recipe. The amount of each component of the recipe is shown in reference to 100 parts of pure monomers. The base rate of each charge stream shows a volume to be charged that is compensated for impurities.

The rate of charging all of the component streams takes as its reference base the time it takes to fill 12 reactor vessels when charging all streams at their base rate.

The thpt. factor (Eq. 43) shows how many times the base rate volumes of all streams that is desired to be charged to the line. The required production determines the thpt. rate charged. Following is the sequence of calculations:

37. Base rate (parts/100 parts monomer)/(purity) (lbs/gal.) 38. Base parts the parts at which the base rate is calculated 39. Target parts desired parts of that individual stream to be charged 40. GPMBI target base rate (base rate/base parts)* (target parts) 41. Volume factor constant 12(3750)/60 (0.95)= 789.474 (cold rubber recipe) 12.5 (3750)/60 (0.95)=822.368 (hot rubber recipe) 42. GPMT (g.p.m. total) volume factor constant/actual throughput Thpt. factor GPMT 2 (43) all char e streams 44. GPMI (g.p.m. for individual stream) (thpt. factor)* (GPMBI) 45. Seconds calibration seconds factor/GMPI 46. Chart settings GPMl/chart set factor In the foregoing equations 26-46:

0,, G G gain factors I(,, K constants M change is modifier parts since last calculation E old error E, new error EOL end of line A specific machine program for carrying out the operations above described on the TRW-330 machine is given in Tables IV and V. As will be obvious from the various tables to follow, octal notation is used in the designation of memory sector cations and in the data in the memory sectors, the latter being scaled to 2 Table IV is a dump of tracks 200, 201, 202 and 205, of the machine which contain the status model up-dating program: the initial instruction for this portion of the program is on track 200, at sector 162 as indicated by an asterisk at the instruction in question. This program blends each variable forward in the reactor chain in real time. Table V is a dump of track 213 and contains the reaction program which calculates adjustments to be made on the basis of the status model. The initial instruction is at memory sector 0, as indicated by the asterisk in the table. Both of these programs use the known operation instruction codes of the TRW-33O machine.

Table VI is a dump from tracks 305 to 307 and 310 to 312, containing typical values for the tracked parameters in each of the reactors 64-1 to 64-21. Table VII is a memory map of these tracks, and Table VIII is a glossary of the abbreviations and terms used in Table VI.

Table IX is a dump from track 43, which accepts the telemetered input coming from the field via the lines 17 and 18" and also contains (together with track 54, Table XII) control set points read out via the lines 20. Table X is a memory map of track 43. As indicated on the table, sectors 024 contain the current temperatures of the reactors as telemetered in from the field. Sectors 40-53 contain set points which are calculated by the program of Tables IV and V and imposed on the temperature controls for the indicated reactors. Sectors 56-64 and 71-77 contain input data telemetered from various instruments in the field as indicated by reference characters, or in the absence of reference characters, manually entered from analytical data via the keyboard 16. The values in these sectors are in terms of centivolts, not yet converted to actual values of flows etc. In order to get values averaged over the 5- minute cycle, the values are read in from the field each minute and added up in these registers to provide values proportional to the average over the 5-minute cycle. The even-numbered sectors from to 136 contain flow control output data for the various ingredients in terms of R values" calculated by the program section from the input data and the desired flow of each of the indicated ingredients (butadiene, styrene, etc.) these R values" for each ingredient being multiplied by the corresponding target partsifor that ingredient appearing in Table XII, track 54 even-numbered sectors -176 and the product applied as a set point upon the flow rate controller for that ingredient through an appropriate one of lines 20'. The remainder of the memory sectors in Tables X and XII are values of the indicated variables either calculated from the input data for intermediate use and/or logging by the program sector, or analytical or base rate requirements entered via the keyboard 16.

TABLE VIII BLENDED VARIABLES FOR EACH REACTOR Bd lbutadiene impurity peak from chromatograph Bd 2butadiene impurity peak from chromatograph Bd 3butadiene impurity peak from chromatograph Bd 4butadiene impurity peak from chromatograph Modifier Activator Soap parts Soap pH .Bd 6-butadiene impurity peak from chromatograph 10. Sty 2styrene impurity peak from chromatograph O percent oxygen Rate-reaction rate (conversions/hour) at 0 F. P.C. Densitypre-charge density R.P.reducing power of activator Sty 3styrene impurity peak from chromatograph Conv-calculated conversion DT-reaction time (hours) I-IDT-reaction time (hours at 0 F.) CATCatalyst 0space reserved for new variables 0space reserved for new variables 0space reserved for new variables 0space reserved for new variables 24. 0space reserved for new variables 25. 0space reserved for new variables 26. 0space reserved for new variables 27. IO CONV-reactor 64-10 conversion 28. HDT*HDT desired when midline latex gets to end of line 29. IOI-IDTreactor 64-10 I-IDT 30. XCAT Ptsincremental catalyst parts added at 92 In the foregoing description, reference has been made in general terms to catalysts, activators, modifiers, etc. In general, the catalysts are free-radical yielding compounds, the activator ingredients are valence-changing metal compounds suitable as redox components, and the modifiers are mercaptan-type compounds, all well known in the art. While the invention has been illustrated in connection with the emulsion copolymerization of butadiene with styrene, it will be understood that the system could be applied to the polymerization of any of a wide variety of monomers, both to form homopolymers or copolymers of one or more of the monomers such as the various conjugated diolefins such as butadiene, isoprene, 2,3-dimethylbutadiene, piperylene and the like; styrene and its homologs such as alpha-methyl styrene, the vinyl toluenes, and the like; acrylic compounds such as acrylonitrile, methacrylonitrile, acrylic and methacrylic acid and esters such as methacrylic acid, ethyl acrylate and the like; and vinyl and vinylidene compounds such as vinyl chloride, vinyl acetate, vinylidene chloride, and the like. Likewise the system is not confined to emulsion polymerization, but may also be applied to suspension and solution polymerization, and using other than free radical catalysts, for instance, anionic and cationic catalysts such as lithium and other metal alkyls, Ziegler catalysts, Friedel-Crafts catalysts and the like. it will also be understood that not all of the controls disclosed herein need be exercised simultaneously; for instance, a very good approximate control may be had with the activator loop alone, although it cooperates in a particularly excellent way with the final temperature trim and catalyst control routines. in such case, the other controls would be left to manual manipulation by the operator, on the basis of his own judgment, aided in some cases by values of the various parameters which can be read out from the computer 14. Likewise various combinations of two or more of the controls-say activator and final temperature trim of the postmidline reactors-can be controlled by the computer and other adjustments left to manual control by the operator.

From the foregoing general discussion and detailed described embodiment of the invention, it will be evident that there has been provided a process and system of control for polymerization systems which will operate with a minimum of skilled supervision, will hold the output properties closely to target, and will make progressive adjustments along the line as disturbances are propagated therealong. The system also makes fullest use of all equipment available, and avoids the cutting out of portions thereof for control purposes.

lclaim:

l. in the operation of a train of stirred reactors wherein polymerizable process material is continuously supplied to the first reactor in the train, each successive reactor is connected to continuously discharge its contents of process material into the next succeeding reactor in the train, polymerized process material is continuously discharged from the last reactor in the train and a programmed digital computer is connected to receive signals corresponding to values of process variables in said reactors,

the method of determining the value of a process variable X in each reactor as that variable changes in that reactor as the result of material balance in the reactor train,

said method consisting in the following steps:

generating and storing into memory sectors of a digital computer digital signals corresponding to numerical values of said process variables in the successive reactors and also storing numerical values of other parameters hereinafter tabulated, said numerical values and corresponding memory sectors being designated as follows, in FORTRAN notation Designation of Memory Sector(s) Valucts) Stored X is a process variable N w a subscript, indicating that the variable stored is the value of that variable in the N one of the reactors in the train NUMBER the total number of reactors in the train DTIME. length of time selected for the interval between updatings of the variables X(N) as a member of an array of subscripted memory sectors one for each reactor, the value of the variable X in the N" reactor in the train the volumetric flow rate in the reactors in the train volume of the N" reactor Storing in the memory of said computer a machine-language program corresponding to the following FORTRAN program D0 10 N=l NUMBER 10 X(N)=X(N)+(F/V)*(X(Nl )-X(N))DTIME Causing said computer to periodically execute said program at intervals of time DTlME, thereby at each execution of the program generating updated digital signal values of the variables formerly stored in the array of memory sectors, and storing said updated signal values back in the array of memory sectors X(N) in place of the digital signals present in the array of memory sectors at the beginning of that execution.

2. Process according to claim 1, wherein the variables X(N) are more particularly concentrations of a component of the polymerizable process material in the several reactors.

3. Process according to claim 1, wherein the values of process variables stored in the array X(N) are more particularly values, in the several reactors, of a thermal history parameter, hereinafter designated l-lDT" of the polymerizable process material such that HDT f-Hdt wherein T temperature t= time of subjection of the polymerizable process material to the temperature T wherein the following further steps are periodically performed measuring the temperature in each reactor, generating a digital signal for each temperature so measured proportional to said temperature, and storing the digital signals together with other digital signals in memory sectors of said digital computer designated as follows, in FORTRAN notation Designation of Memory Sector(s) Value(s) Stored storing in the memory of said computer, in addition to the program set forth in claim 1, a further machine program corresponding to the FORTRAN program causing said computer to periodically execute the above program at intervals of time DTlME, in addition to executing the program of claim 1 to update the values of X(N) for materials balance, thereby at each execution of the programs generating digital signal values of X(N) updated additionally for genera tion in the reactors, and storing said additionally updated signals back in the array of memory sectors X(N) in place of the digital signals stored in the array of sectors at the beginning of that execution.

4. Process according to claim 3, wherein the rate of flow of activator to the first reactor in the train is controlled in response to digital signals generated and stored in accordance with the process of claim 3, the further steps which consist in conducting upon these digital signals the summation thereafter subjecting these digital signals to the manipulation Tncw old thereafter subjecting these digital signals to the manipulation reaclurML ML ML/ I reactor:

thereafter further subjecting these digital signals to the manipulation Hdz HDT(N1)= value of the integral f l-ldt for the contents of the (N1)th reactor T(N) temperature in the Nth reactor T temperature in degrees Fahrenheit t= time of subjection of reactants to temperature T Equations (II) and (III) are in terms of FORTRAN N as a FORTRAN subscript, ascribes the variable which it subscripts to the contents of the Nth reactor in the train of reactors n as a subscript, ascribes the variable which it subscripts to the contents of the nth reactor in the train of reactors C the projected final conversion of the monomers charged R the target basic rate of conversion, referred to 0 F., the sub-subscripts new" and old indicating the value of this target rate before and after the operation indicated by equation (2) H value of H in the nth reactor 1 holdup time in the reactor in question X has the value (70C,,/20)

C, conversion in the nth reactor r has the value indicated by equation (2) and is the apparent reaction rate up to midline, referenced to 0 F.

C is the actual conversion, measured at midline, of the reactor train by means of a densitometer In equation (6), the integral expression was evaluated by computations of claim 7, integrated up to the midline reactor A value defined by equation (7) applied as a control to the flow of activator No. l of 16 PMs-5c U 'ITED STATES PATENT OFFICE (ED/EL" CERTIFICATE 0F CORRECTION Pa tent No 3 61 4- 682 Dated 00 tober l9 19?].

Inventor) Robert C. Smith It is rertified that error appears in the above-identified patent and that 531d Letters Patent are hereby corrected as shown below:

ca. 8 11a. 53 90101-31, 2 0 lhmlld "Id "DO 10 1-135-- lhould road 12) 1m (10 m (I) Mum 01-1) m m)- 1 0 I8 (I) I W 11.11. 72. "1th" should and --tht-- :01. 1 11:10 39, 19.0 18mm) 20)' :hculd rad -(9.o um): 20)

11a. 66, 9.oucr(mn)*20)" mum rand (9.0 Acrrxzcnhzo) 11m 70, the 00nd "I" mould rand n 4 ca. 1 line 5 anon". both occurrences, mum read --(:nx.

Q01- 15, 11m 01, the ('03) Z: it the right inltufl 01' to th. 10!:-

11110 +5, "GMPI" nhould rnld --GPllI-- 11m 9+, "EOE. 0nd of 11:10" lhuuld rnfl IBI- I and of 11:10-- line 60, "nailed to 2 should rad --lcllid to 2' 001. 16, .hbln IV, V, VI Ind VII have been omittod Table: II, I Ind XXX have bun omitted (tar 11a 70) C01. 18, 11a. 62 "11 I!!!) I I0!) 'm'n-m" would and .41 x0?) I 10!) MI) mum-- PO-IOSO (5/81)- Patent No.

UNITED STATES PAT FNT OFFICE Dated Inventofls) It: is certified that: error appears in the aboveand that said Letters Patent are hereby corrected as Robert C. Smith October 19 1971 Column 16 after line 3 insert 0027-000 0200-053 0017-001 0000-001 0000-001 0200-051 0027-036 0117-001 0025-000 0075-OQO 0117-001 0200-136 0117-001 0117-001 0117-001 0200-065 0056 0117-901 0056 0117-001 00%- O2Sk-0O7 0056 0117-001 0117-001 0200-052 0021-036 0117-001 0017-001 0117-001 0117-001 0027-000 020 0-130 0200-157 TABLE IV BIAI'UB mm. PROGRAM DUMP 0200-053 0077-000 OO2-v-00O OOZHOOO 0060-0 51 OOZJ-OOU 0117-001 0201-101 0055-023 0117-001 0201-01 0 0021-000 02510-002 OOZh-DMO 0201-051 0117-001 OOZh-OOG 0070-017 2650-020 0665-076 0117-001 0200-107 0052-077 0117-001 0202-17 0052-131 GOW 0102 0076 0076 0056 1025 identified patent shown below:

0030-000 0030-000 (28-06-160 0200-(5'7 0117-001 OCLZk-OOD 0025-000 0200-123 0000-177 0025-172 0117-001 0117-001 C117-00l 33 35 No. 3 of 16 UNITED STATES PATENT OFFICE CERTIFICATE OF CORRECTION Patent No. 3,615,682 Dated October 19. 1971 Inve t t C It is certified that: error appears in the above-identified patent and that: said Letters Patent are hereby corrected as shown below:

F- Column 16, after line 3 and after Table 1V(lst part) ,insert .1

ILBL! IV (CONT'D.)

IDQOZ-WO 0056 0117-001 01'02 0029-000 0056 0117-001 0030 0000-001 0056 0117-001 02 42 0025-000 0056 0117-001 1071 0062-06 0000 0200 0000 0201-177 0000 0202-1 6 0000 0201-077 0056 0117-001 0000 0000-001 0000 0000-002 0000 coon-00 0000 0000-010 0000 0000-020 0000 0000-0 0 0000 0000-100 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 C056 117-001 0056 0117-001 0056 011 -001 00% 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0055 0117-001 0056 0117-001 0056 o117-oc1 0056 0117-001 0056 0117-0-21 0056 0117-001 0056 0117-001 0052 0006-000 0070 0052-022 0-20 0001-11 7 00% 0225-057 0056 0171-002 0000 0225-071 02- 0000-025 001k 0225-063 ooh: 0030-000 0000 0225-071 00 0 0075-000 0 0000-001 @050 0115-062 0076 0052-022 0076 0005-167 0000 0225- 00 2 0006-000 001k o222-o 2 10 2 0006-010 0000 0225-075 0056 0117-001 ooco 0225477 0055 0117-001 0006 0215-152 01 2 oozh-ooo 0056 0117-001 0076 006 -10 0-20 0001-177 nous 0063-112 0056 0171-002 0070 0077-000 0-2 m o-0 5 0033 0063-133 00-2 0030-000 0 0 0000-012 0030 0075-000 02-2 0025-000 0050 0115-062 006 033M177 00 s Dom-167 0000 0225- 0061 0252-1 7; 0036 0052-0-- o2 +2 0025-000 0000 0227-071 0071 033M177 0056 0117-001 ony's 0065-16? 0 .26 0000-002 0-70 0115-062 owl 05-477 0 32 0225-1 -1 00%;: 0035-000 0056 0070-000 0050 0222-137 0100 0110-000 0000 0213- 1200 02254-3 00% 0117-001 0024 0365-167 0056 0117-001 0-70 02-0-07 00% 0117-001 mm: 0110-012 0056 0117-001 0000 ooao-om 0056 0117-001 005- 0215-060 0056 0117-001 0000 0211-000 0055 0117-001 0056 01:7-c01 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0050 0117-001 005:: 0117-001 00% 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 0056 0117-001 005:) 0117-001 Po-mso (5/69;-

Patent No.

Inventor(s) UNITED STATES PATENT OFFICE No of 16 Dated Robert C. Smith October 19. 1 9' 7L It is certified that error appears in the above-identified patent and that said Letters Patent are hereby corrected as shown below:

Column 16 after 1 ine 3 and after Table IV cont d insert TABLE V "ACTION PEUGRM DUMP 0000-001 0000-022 O11" 001 CC? -05'\ 0017-001 0095-000 0213-005 0070-000 0117-001 0252-105 0117-001 011 T-OO]. 0117-001 0025-000 1352-161 0117-001 N0. 5 of 16 232 3 3 UNITED STATES PATENT OFFICE CERTIFICATE OF CQRRECTIQN Patent No. 3 ,61+,682 Dated October 19 1 1971 Robert C. Smith Inventor(s) It. is certified that error appears in the above-identified patent and that said Letters Patent are hereby corrected as shown below:

"column 16, after 11116 3 and after Table v, insert IABLE VI TYPICAL DUMP OF TRACKS 305-? AND 310-312 0000365613 0002153 55 00003 2100 0000176767 0003350317 0000000001 0002122503 0005200006 0001127007 0001251000 00052 3650 0000500575 000 0000 000 365207 0002553 25 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000365607 0002153555 00003 2100 0000176762 000 62h310 0000300001 0002121557 0005202733 00011270 0 00012530 5 00046 5306 0000570575 000 00000 0003652030 00025556% 000153222 0000165671 0001160027 0017736703 1777777773 1777777773 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000365602 0002153655 0000332100 0000176755 000562 +260 0000000000 0002121310 0005206267 0001127007 0001257130 0004303175 0000550575 0000000000 0003651765 0002560156 0003031266 000032h33k 0002101053 0017735656 1777777766 1777777766 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0002153 65 0000332100 0000176750 000 621161 0000500000 0002121072 OOOSZlOlkfi 0001127067 0001261 427 000%26572 00005 0575 000 0000 0003651720 0002560555 000 1311723 0000455022 00031505 6 0017765621 1777777761 1777777761 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 000371 06 31 0000531k63 0003665107 0000365571 0002153 655 0000332100 0000176703 000 620025 0000000000 0002120750 000521027 0001127057 0001261657 0006176053 00005 4766 00 0000 0003651663 0002557555 00056 0273 0000623521 000020 31 0020020557 1777777757 1777777757 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000365571 0002153055 0000332100 0000176736 0005623605 0000000000 0002120656 0005210311 0001127057 0001261320 0000156650 0000562766 OOOMOOOOO 0003651617 0002556150 00077175 03 0000763251 0005225013 0020060515 1777777757 17777777'+7 00 0000600 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 00000000c0 0000000c00 0000365 5 71 0002153060 0000 332100 0000176731 000562 33 5 0000000000 0002120666 0005210153 0001127053 0001260361 00051 1620} 000094766 000 2100000 0003651906 0002555055 0010553722 0001122766 000627 6070 0020112023 l7777777 +2 17777777 2 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0O0O*Ol7? 0002271633 00003 +6216 0000163117 0005671051 0000000000 0002131277 0005232217 0001117036 00013306 7 0002556520 0000502766 0000000000 0003636507 0002671613 0013220232 0001352113 0007751513 0021150621 1777777731 1777777731 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 N0. 6 of 16 Egg? UNITED STATES PA'IIINT OFFICE .CERTIFICATE OF CORRECTION Patent No. 3,610,682 Dated October l9. 1971 I r( Robert C. Smith It is certified that erro r appears in the above-identified patent and that said Letters atent are hereby corrected as shown below:

Column 16, after line 3 and after Table VI, installment beginning "A0 305-000" insert TABLE VI (CORTD) 0000365571 0002153671 0000 2150 0000176720 000 623217 0000000000 0002120 12 00052101 2 0001127323 0001257750 000012 00 0000502 0 000 ;700:100 0003651 735 0002556357 00130162 1 0001022263 0010337262 00201 3 3222 1777777731 1777777731 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000 65612 0002150610 0000 32025 00001767 0 000562 07 0000000000 00021207 0005210 00 0001130020 00012607 0 00000532 1 0000502 0 000 0300000 000 65131 00O25627 0 O0l62 r756l 0001721575 0012377235 00201276 1 1777777720 1777777720 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000365632 000215535 0000 26 0 O00Ol767 000 623%75 0000000000 0002120753 0005210022 00011 1270 00012616 1 0000022 7 0000502 000 0000 0003651253 0002505720 0017727025 000 1 12 0013 77373 0020120521 1777777713 1777777713 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0000000000 0002156 0 0000 126 000017677 000562 605 0000000000 0002120752 0005210555 00011 2517 0001262 76 0003769075 0000500556 000 000000 000 651172 0002571 02 0021 72122 00022210 0 001M560071 0020105272 1777777706 1777777706 0000000000 0000000000 0000000000 0000000000 0000000000 0017669130 00322 065 0013376320 0000000000 0000000000 0000000000 0000000000 10310-000 0000362337 0002157507 0000333537 0000177055 000 623776 0000 00000 0002120736 0005210711 71 00012903 1 0003725510 0000537015 000 000 0003651121 0002575322 002303 4567 00023605 5 00157155 0020067101 1777777701 1777777701 0000000000 0000000000 0000000000 0000000000 0000000000 0017627 315 00 32 3230-61 001 3370723 0000000000 0000000000 0000000000 0000000000 0000366M 0002161167 0000330301 0000177202 000569 055 0000400000 0002120711 0003511025 0001136120 0001266011 0003660111 0000536066 0000000000 0003651050 0002001532 002 775 1 0002520260 0017051007 0020056677 1777777670 177 777767 0000000000 0000000000 0 0000000000 0000000000 001757665 00 32 3721 75 001 3373200 0000000000 0000000000 0000366172 0002163122 0000335205 006017 7 0005623772 0000 200000 0002120651 00052111 12 00011 055 7 0001267561 0003622000 0000535705 0000000000 0003650777 0002606073 0036100510. 000266000 0020205353 0020025122 1777777667 1777777607 0000000000 0000000000 0000000000 0000000000 0000000000 0017553010 00320 6677 0013 7l 56 0000206336 0000000000 0000 66 60 0002165300 000033629 0000177700 000 6239A 0000 0000 0002120600 0005211226 00011 73600 0001271370 0003562060 0000535215 00 00000 0003650726 0002612572 0027612030 0003017530 0021301370 00200021 3 1777777662 1777777662 0000000000 0000000000 0000000000 0000000000 0000000000 00175322 0032072725 0013370002 0000502021 0000000000 0000000000 0000000000 

2. Process according to claim 1, wherein the variables X(N) are more particularly concentrations of a component of the polymerizable process material in the several reactors.
 3. Process according to claim 1, wherein thE values of process variables stored in the array X(N) are more particularly values, in the several reactors, of a thermal history parameter, hereinafter designated ''''HDT'''' of the polymerizable process material such that HDT Hdt wherein H 2T/18 T temperature t time of subjection of the polymerizable process material to the temperature T wherein the following further steps are periodically performed measuring the temperature in each reactor, generating a digital signal for each temperature so measured proportional to said temperature, and storing the digital signals together with other digital signals in memory sectors of said digital computer designated as follows, in FORTRAN notation Designation of Memory Sector(s) Value(s) Stored X(N) as a specific instance of X(N) referred to in claim 1, and as a member of an array of subscripted memory sectors, one for each reactor, the value of HDT for the contents of the Nth reactor in the train T(N) as a member of an array of subscripted memory sectors, one for each reactor, the value of the temperature in the Nth reactor in the train storing in the memory of said computer, in addition to the program set forth in claim 1, a further machine program corresponding to the FORTRAN program DO 11 N 1, NUMBER H(N) 2 ** (T(N)/18) 11 X(N) X(N) *DTIME causing said computer to periodically execute the above program at intervals of time DTIME, in addition to executing the program of claim 1 to update the values of X(N) for materials balance, thereby at each execution of the programs generating digital signal values of X(N) updated additionally for generation in the reactors, and storing said additionally updated signals back in the array of memory sectors X(N) in place of the digital signals stored in the array of sectors at the beginning of that execution.
 4. Process according to claim 3, wherein the rate of flow of activator to the first reactor in the train is controlled in response to digital signals generated and stored in accordance with the process of claim 3, the further steps which consist in conducting upon these digital signals the summation thereafter subjecting these digital signals to the manipulation
 5. RT RR *60/Cs thereafter subjecting these digital signals to the manipulation thereafter further subjecting these digital signals to the manipulation
 7. A0 AML*(RT/rML) and adjusting the flow of activator proportionally to the finally generated signal Ao, the notation of the mathematical expressions hereinabove being as follows: DTIME interval of time for each cyclic repetition of the calculations HDT the value of the integral as specified in claim 3 NUMBER number of reactors in the train of reactors HDT(N) value of the integral Hdt for the contents of the nth reactor F rate of flow of reactants V reactor volume HDT(N-1) value of the integral Hdt for the contents of the (N-1)th reactor T(N) temperature in the Nth reactor H 2T/18 T temperature in degrees Fahrenheit t time of subjection of reactants to temperature T Equations (II) and (III) are in terms of FORTRAN N as a FORTRAN subscript, ascribes the variable which it subscripts to the contents of the Nth reactor in the train of reactors n as a subscript, ascribes the variable which it sUbscripts to the contents of the nth reactor in the train of reactors Cs the projected final conversion of the monomers charged RT the target basic rate of conversion, referred to 0* F., the sub-subscripts ''''new'''' and ''''old'''' indicating the value of this target rate before and after the operation indicated by equation (2) Hn value of H in the nth reactor Tau holdup time in the reactor in question Xn has the value (70-Cn/20) Cn conversion in the nth reactor rML has the value indicated by equation (2) and is the apparent reaction rate up to midline, referenced to 0* F. CML is the actual conversion, measured at midline, of the reactor train by means of a densitometer In equation (6), the integral expression was evaluated by computations of claim 7, integrated up to the midline reactor Ao value defined by equation (7) applied as a control to the flow of activator 